Effects of the Sound Speed of Quintessence on the Microwave Background and Large 

Scale Structure 
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We consider how quintessence models in which the sound speed differs from the speed of light 
and varies with time affect the cosmic microwave background and the fluctuation power spectrum. 
Significant modifications occur on length scales related to the Hubble radius during epochs in which 
the sound speed is near zero and the quintessence contributes a non-negligible fraction of the total 
energy density. For the microwave background, we find that the usual enhancement of the lowest 
multipole moments by the integrated Sachs- Wolfe effect can be modified, resulting in suppression 
or bumps instead. Also, the sound speed can produce oscillations and other effects at wavenumbers 
k > fO~^ h/Mpc in the fluctuation power spectrum. 

PACS numbers: PACS numbor(s): 98.80.-k, 98.70.Vc, 98.65.-r 



One of the greatest challenges in cosmology today is 
to identify the nature of the dark energy component 
that comprises most of the energy density of the uni- 
verse and that is causing the expansion of the universe to 
accelerate. jl| Two candidates are a cosmological constant 
(or vacuum density) and quintessence, a dynamical en- 
ergy component with negative pressure. Distinguishing 
the two is important for cosmology in order to refine our 
knowledge of the composition of the universe and to trace 
more accurately its evolution. It is even more important 
for fundamental physics since it informs us how we must 
modify unified theories to incorporate dark energy. 

One way to distinguish whether the dark energy is due 
to a cosmological constant or quintessence is to measure 
the equation of state, w, the ratio of the pressure p to 
the energy density p. A cosmological constant always has 
w — ~l whereas a scalar field generally has a wiz) that 
differs from unity and varies with red shift z. Through 
measurements of supernovae, large-scale structure and 
the cosmic microwave background (cmb) anisotropy, the 
equation of state may be determined accurately enough 
in the next few years to find out whether w is different 
from —1 or not. 

A second way to distinguish the nature of dark en- 
ergy is to measure its sound speed to determine if it is 
different from unity. The sound speed can be detected 
because it affects the perturbations in the quintessence 
energy distribution. This approach is less generic be- 
cause the sound speed in many models of quintessence in 
the literature is equal to unity (the speed of hght), e.g., 
models in which quintessence consists of a scalar field (0) 
with canonical kinetic energy density {X = \ {d^(f)f') and 
a positive potential energy density {y{4>)). However, in 
general, the sound speed can differ from unity and vary 
with time. Detecting these effects is an independent way 
of showing that dark energy does not consist of a cosmo- 
logical constant. 



An important motivating example is fc-essence.|3j In 
these models, the A;-essence undergoes two transitions 
in its behavior, one beginning at the onset of matter- 
domination and a second when fc-essence overtakes the 
matter density. During the radiation-dominated era, the 
fc-essence energy tracks the radiation, falling as 
where a is the scale factor. This tracking feature is sig- 
nificant because it explains how dark energy can remain 
a subdominant but not completely negligible component 
prior to matter-radiation equality. Generic models with 
this behavior (fc-essence is only one example) are known 
as "trackers." |^ The tracking property is relevant to our 
purpose here because the sound speed can have greater 
influence if the dark energy is non-negligible for a broad 
range of red shift. (A cosmological constant becomes 
completely negligible by red shift z = 2, but this is not 
the case for typical tracker-type models.) 

What is special about A;-essence models is that the on- 
set of the matter-dominated era automatically triggers a 
change in the behavior of fc-essence such that it begins to 
act as an energy component with w{z) w —1. The effect 
is achieved with a scalar field by introducing a kinetic en- 
ergy density which contains higher order derivative terms 
(that is, terms non-linear in X) and which results in at- 
tractor solutions with the desired behavior. This kind 
of model is appealing because it explains dynamically 
why cosmic acceleration only begins well after matter- 
domination, just as we observe, which is one of the puz- 
zling aspects of dark energy. When A;-essence overtakes 
the matter density, w(z) changes to another value greater 
than —1, the precise value of which depends on the de- 
tailed model. The way in which the attractor behavior of 
fc-essence is automatically controlled by the background 
equation of state is important motivation for cosmolo- 
gists to consider /c-essence models, specifically. However, 
the analysis in this paper is more general. 

For the purposes of this paper, what is significant 
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FIG. 1: The speed of sound as a function of the scale factor 
1 (ffltoday = 1) for a A;-essence model studied in Ref.0]. Note 
the changes in Cs between the last scattering surface and the 
present epoch triggered by the transformations in the back- 
ground equation of state. 



about the /c-essence example is an indirect consequence 
of the tracker behaviour: the sound speed Cg undergoes 
dramatic changes as w{z) passes through a series of tran- 
sitions in evolving from the radiation- to the matter- to 
the fc-essence-dominationed epochs. The behaviour of 
Cs in a particular fc-essence modelQ is shown in Fig. ^ 
Here the sound speed (i) is relativistic in the radiation- 
dominated epoch, (zi) has rapid changes between zero 
and relativistic as ■w{z) settles towards the attractor so- 
lution, which corresponds to c| ^ and w{z) — > —1, 
and, then, {in) increases somewhat at recent red shift 
when /c-essence overtakes the matter density. The equa- 
tion of state in tracker models must always have these 
three phases. Without knowing what the dark energy 
consists of, it is reasonable to imagine a generic class of 
dark energy models where the sound speed during each 
phase varies in a manner similar to the specific model 
described here. In particular, changes in sound speed 
are related to the changes in the scalar field as it moves 
rapidly between the radiation tracking and attractor so- 
lution in phase (ii), and asymptotes to the w = —1 at- 
tractor point {Hi). 

Our study here does not rely specifically on non-linear 
kinetic energy or other aspects specific to fc-essence, but 
simply on the equation of state w{z) and the sound speed 
Cs{z) as a function of red shift z. The purpose of this pa- 
per is to study possible efl^ects of the sound speed and 
its time-variation on the cosmic microwave background 
anisotropy and on large scale structure. In an earlier 
paper, the subtle changes due to sound speed in the 
position and shape of the acoustic peaks in the CMB 
power spectrum were investigated and it was shown that 
the effect of the sound speed could be discriminated from 
other physical effects provided that the fc-essence den- 
sity is non-negligible (greater than 1%) at last scatter- 
ing. Here we focus on the large angular scales corre- 



sponding to the quadrupole and first few higher multipole 
moments. We show that ^ 1 at small red shift can 
produce a suppression of these first few multipole mo- 
ments, an effect similar to the low quadrupole observed 
in the COBE data. 0, Q We also show how rapid vari- 
ations in the sound speed, as in the A:-essence example, 
can produce bumps, oscillations and other novel features 
at small scales in the mass power spectrum. The precise 
effect depends on the detailed model. 

As noted above, our calculations and examples do not 
assume that the dark energy consists of fc-essence specif- 
ically, which is a purely kinetic action for the scalar 
field. In fact, our computer code, a modification of the 
standard CMBFAST program, 0, ^ effectively treats the 
quintessence component the same as a scalar field with a 
potential. The code automatically determines a potential 
that produces the desired w{z). Independently, the per- 
turbation equations are modified to include the desired 
sound speed, Cs{z). 

The effect of the speed of sound on the CMB perturba- 
tion equations is such that for <C 1, fc-essence fluctu- 
ations are enhanced via gravitational instability by the 
cold dark matter (cdm) potentials. For those familiar 
with the code, the modifications are straightforward. The 
perturbed line element is 



(1) 



where Sij is the unperturbed spatial metric, and hij is 
the metric perturbation. We shall use h to represent the 
trace of the spatial metric perturbation. The effect we 
are examining is due to the perturbations to the fc-essence 
stress-energy in the synchronous gauge for a mode with 
wavenumber fc (we omit the fc subscripts here): 

dp = -2p— - {p + p)—Cs^ 

<p y 



Sp 



-2p 
1 



d(t) 



(P + P) 



5y 



(2) 



where y= l/\/X and 9 is the divergence of the fiuid ve- 
locity. The density contrast, 5 = bpjp, obeys the equa- 
tion 



-(1 + w) 



1 



6p 

Tp 



w S, (3) 



where the derivative is with respect to conformal time 
and Ti. = d/a. The relation between Sp and Sp, derived 
using equations ^ above combined with the background 
equation of motion, is 



Sp = clSp + -| [3n{l + w){cl -w)+w], 



(4) 



which leads to a simplified evolution equation for the 
velocity gradient 

9 = {3cl~i)ne + cies/{i + w). (5) 
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FIG. 2: Comparison of the lowest multipole moments of the 
CMB temperature power spectrum for a series of models with 
w = —0.8: (a) Cs = 1 (dotted); (b) Cs = 1 until 2 = 5 and then 
Cs = for z < 5 (solid); and (c) Cs = for all z (dot-dashed). 

We can see in equation (|3J| that a small sound speed will 
cause the velocity gradient to decay; with the conven- 
tional gauge choice, 9cdm = 0, the inhomogeneities in the 
fc-essence will describe a fluid which is comoving with 
the cold dark matter. From equation Q), we see that the 
second term on the RHS will be negligible even on scales 
approaching the horizon. The overall effect is that the 
pressure fluctuations 5p are weak and fc-essence fluctu- 
ations are enhanced via gravitational instability by the 
CDM gravitational potentials. The CMBFAST code takes 
w{a) and Cs{a) as inputs, so it is possible to manually 
adjust these functions to have any values (including, of 
course, Cs = 1). 

As a first set of cases we consider a sequence of models 
in which the equation of state is constant, w — —0.8, 
and we vary the sound speed. Note that this example 
does not correspond to tracker-type evolution and, like 
the cosmological constant, the quintessence component 
becomes negligible at 2; > 2 or so. Fig. [21 compares three 
cases: (a) = 1 (dotted); (b) = 1 for z > 5 and then 
Cs = for z < 5; and (c) = for all z. All models have 
Vlrn = 0.25, Hq ^ 70 km/s/Mpc, and Qb = 0.04, where 
^m(6) is the ratio of the matter (baryon) density to the 
critical density and Hq is the Hubble parameter. 

We observe the following features: 
(1) The (late-time) integrated Sachs- Wolfe (ISW) effect 
is suppressed if the sound speed is near zero at recent 
red shift. The ISW effect occurs in models with < 1 
because the gravitational potential decreases as photons 
propagate between the last scattering surface and to- 
day. The net effect is an increase in the multipoles on 



angular scales which enter the horizon when flm < 1, 
that is to say, the \ow-£ multipole moments. This is 
the effect predicted for models with a cosmological con- 
stant, for example. Dark energy with sound speed near 
zero suppresses this effect because the fluctuations in the 
quintessence component add to the gravitational poten- 
tial and compensate for the decrease one would obtain if 
the quintessence were smoothly distributed. jSl] 
(2) For constant w models the ISW suppression is only 
weakly dependent on the time-dependence of the sound- 
speed. For example, the lower two curves in Fig. [21 
nearly overlap. The reason for this is apparent - the 
quintessence component is a negligible fraction of the en- 
ergy density for z > 2, so the behavior of the sound speed 
is only important over a narrow range of recent red shift. 

We have not shown a plot of the multipole moments for 
large values of £ because the curves are essentially degen- 
erate over this range. Similarly, the mass power spectra 
for these models are virtually identical. Hence, the only 
effect of the sound speed for these non-tracker type mod- 
els with nearly constant w is a modest suppression of the 
low-£ multipole moments in the CMB. 

In Fig. [3 we show the small-i? CMB multipole mo- 
ments for models with the same u;(z), where w is now 
z-dependent. As a specific example, we choose w(z) 
in these examples to correspond to the fc-essence model 
shown in Fig. ^ and discussed in Ref.0]. As before, 
all models have current values of ^Im = 0.25, Hq = 
70 km/s/Mpc, and fif, = 0.04. The details of w{z) are 
not important here except that it has the feature that the 
dark energy is a tracker component that is non-negligible 
at last scattering and earlier, roughly 10% of the critical 
density. We note the following features: 

(1) Comparing the models with Cg = 1 (dotted) and 
with the case with = 1 for z > 10 but Cg = for 
z < 10 (solid), there is greater suppression than in Fig. [21 
because the quintessence component is non-negligible for 
a broader range of red shift, enhancing its clustering ef- 
fect. (We have considered this example of time- variation 
because it is generally plausible that the sound speed of 
the quintessence component changes as it overtakes the 
matter density, as occurs in the case of fc-essence.) 

(2) More general time-variation of the sound speed 
can produce bumps and wiggles in the \ow-£ multipole 
moments. The remaining examples in Fig. [3 are cases 
where the sound speed corresponds to the actual fc- 
essence model shown in Fig. ^ (dot-dashed), or simple 
variations. For the variations, we have considered cases 
where the sound speed is = 1 until last scattering, then 
drops to Cs = except for a spike near z = 100 where 
the quintessence component approaches the attractor so- 
lution with w = ^1. While this roughly mimics features 
that are characteristic of fc-essence, here we have bounded 
Cs < 1 and treated the spike as a smoothly varying func- 
tion to illustrate that the unusual features seen in the ac- 
tual fc-essence example do not depend on having a sound 
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FIG. 3: Comparison of the lowest multipole moments of the 
CMB temperature power spectrum for a series of models with 
the same w{z) (in this case, corresponding to the fc-essence 
model in Fig. and Ref.^) but different Cs{z): (a) Cs = 1 
(dotted); (b) = 1 for z > 10 and Cs = for z < 10 (solid); 
the fc-essence model shown in Fig. Q (dot-dashed) ; and two 
variations (short- and long-dashed) with rapid variations in 
(spikes) as described in the text. 



speed greater than unity or with near-discontinuous be- 
havior. We find that the precise form of the low-£ power 
spectrum is sensitive to the red shift and width of the 
spike, producing in some cases bumps and wiggles in the 
small-i? multipole moments. 

Curiously, the COBE data suggests a suppression, 
and perhaps even a bump-like feature in the low-£ 
multipoles.]^, The data is not decisive because the 
cosmic variance and experimental error are large for the 
\ow-£ multipoles. Forthcoming data from the MAP and 
Planck satellite experiments may lead to some improve- 
ment. 

Fig.^shows the full CMB power spectra for the same 
set of models as in Fig. O Recall that these models 
are completely identical except for differences in c^(z). 
These spectra are roughly similar, but could be easily 
distinguished in near future CMB power spectrum mea- 
surements, e.g., by the MAP and Planck satellites. How- 
ever, detecting and distinguishing c'^{z) is not as sim- 
ple as the Fig. 0] may suggest. The modest differences 
in peak heights and positions seen here can be approx- 
imately mimicked by varying other parameters. So, the 
challenge becomes to separate effects of (z) from those 
of the other parameters. As shown in Ref.[2|, the effects 
are not completely degenerate, so highly precise measure- 
ments of the power spectra can make the separation, in 
principle. However, the differences are small enough that 



a systematic search through possible forms for (z) along 
with the standard parameters is required in most cases 
to detect them, and this is not part of any current fitting 
procedure. Observing anomalies in the low-£ multipole 
moments and in the mass power spectrum (see below) as 
discussed in this paper would provide strong motivation 
for developing such a procedure. 

We note that similar suppression of the \ow-£ multi- 
poles has been obtained previously by keeping = 1 
but varying the equation of state w(z).01 The greater 
is the value of w, the greater are the fluctuations in the 
quintessence field, which is the same trend one obtains 
by lowering the sound speed. However, if = 1, the 
suppression of multipoles only occurs if w{z) is greater 
than —0.2 today. This is inconsistent with the ob- 
served acceleration of the universe, which requires w{z) < 
— 1/3.0,01 On the other hand, keeping w(z) « —1 today 
and varying the sound speed history is an interesting al- 
ternative because it is physically motivated and does not 
run into conflict with any other current observations. 

Similarly, one could modify the low-^ multipoles by 
altering the primordial power spectrum, for example, 
by introducing an inflation potential with bumps or 
dips strategically placed so as to produce suppression or 
bumps in the low-^ multipoles. This solution is unappeal- 
ing. Introducing such features at any point in the inflaton 
potential requires very awkward fine-tuning; but, in this 
case, where we are seeking to alter only the low-i multi- 
poles, such an approach is even more awkward because 
the features in the inflaton potential must be introduced 
in a way so as to affect specifically those wavelengths 
that recently entered the horizon (rather than longer or 
shorter wavelengths) . There is simply no good reason for 
this. In contrast, it is natural to associate the suppression 
or bumps in the low-£ multipoles with dark energy be- 
cause it only dominates at late times and so will naturally 
affect modes whose wavelengths only recently entered the 
horizon. 

Fig. El compares the total ffuctuation power spectrum 
(including contributions from ffuctuating components) 
for the same models as in Fig. |3 We have intention- 
ally included an example with exaggerated features to 
illustrate that the sound speed effect can be quite large 
(in this case, by making an optimal choice of the red 
shift and width of the spike in Cs). The main effect to 
be observed here is that there can be enhancement in 
power, suppression on some scales, oscillations and other 
features in P{k) at large k > 0.01 h/Mpc. Clearly, un- 
certainty in the sound speed can lead to non-negligible 
uncertainty in erg and, more generally, the shape of the 
power spectrum even when all other parameters are fixed. 
The effect could be important in explaining anomalies 
that may arise in P(fc) as high quality red shift survey 
data is obtained and compared to CMB measurements. 

When = 0, the quintessence fluctuations are en- 
hanced, producing the effects in the microwave back- 
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FIG. 4: Comparision of higher multipole moments of the 
CMB temperature power spectrum for the models in Fig. |3 
The spectra have been normalized so that the amplitudes 
match at the top of the first acoustic peak. 
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tion is modified, 

^dH 3,-2 ,d6 feci 2 A . n 

where a is the Friedmann-Robertson- Walker scale fac- 
tor and A[cl,'w\ = 1 — 5w + 2c^ and B[cl,w\ — 
1 - 6c2 + 8w - Sw^).!!!! For = w = 0, we have 
A[c^,w;] = B[(?g,w\ = 1 and we obtain the conven- 
tional dust-dominated growing solution, 5 a. How- 
ever, in general, the "growing" solution is 5 oc a''' where 
7 = i (1 - I A -h [(1 - -f 6B] 1/2) ^ which is less than 
zero for w < —0.12 (for all > 0). Hence, despite the 
fact that — 0, the accelerated expansion dominates 
and causes the fluctuations to decay. In a universe with 
a mixture of matter and quintessence, the matter slows 
the expansion, so the quintessence may cluster when it 
is subdominant. Indeed, it is the gravitational influence 
of matter fluctuations on quintessence that is enhancing 
the perturbations in quintessence that have been studied 
here. However, this clustering is halted once quintessence 
begins to dominate and the effective value of w becomes 
negative. The only significant effects are the kinds of 
changes in the microwave background and mass power 
spectra seen here. 

What we conclude, then, is that modest changes in 
Cs{z), as motivated by some existing models of dark en- 
ergy, can produce modest but measurable changes in the 
CMB and mass power spectra. The highly precise data 
obtained from the MAP and Planck satellites and from 
the Sloan Digital Sky Survey may reveal these subtle ef- 
fects. The precise behavior of the sound speed is, by 
itself, of limited interest. But, what is important is that 
the detection of any deviation from = 1 would be a 
direct sign that the dark energy is a complex, dynamical 
ffuid rather than an inert cosmological constant. Hence, 
it is a target well worth pursuing. 
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by and NSF Graduate Research Fellowship (SD), NSF 
grant PHY-0099543 (RRC) and Department of Energy 
grant DE-FG02-91ER40671 (PJS). 



FIG. 5: Comparison of the shape of the total fluctuation 
power spectrum, P{k) as a function of wavenumber k for the 
sequence of models in Fig.|21 The normalization of the curves 
is arbitrary. 



ground and mass power spectrum that we have observed, 
but the ffuctuations do not cluster like matter because 
the equation of state, w, is different. For example, con- 
sider a universe containing only quintessence with = 
and perturbations S = Sp/ p, where p is the quintessence 
energy density. Since w ^ 0, the Jean's instability equa- 
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